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Abstract 

We study noise spectra of currents through a tunnel junction in weak tun- 
neling limit. We introduce effective capacitance to take into account the 
interaction effect and explicitly incorporate the electromagnetic environment 
of the junction into the formulation. We study the effect of charging energy 
and macroscopic environment on noise spectra. We calculate current fluctu- 
ations at tunneling barrier and fluctuations measured at leads. It is shown 
that two fluctuations have different noise spectra and the relation between 
them is nontrivial. We provide an explanation for the origin of the difference. 

Experimental implications are discussed. 
PACS numbers: 72.70.+m, 73.40.Gk, 73.40.Jn 
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I. INTRODUCTION 



Even though "noise" usually represents unwanted fluctuations which blur information, 
noise is on the fundamental level closely connected with the dynamics of a system and it 
contains variety of information. 

In this context, there have been many recent studies on a noise spectrum of current 
in small devices. In ballistic systems, the noise power can be suppressed to zero due to 
Pauli exclusion principle In diffusive systems, the noise power is suppressed to 1/3 

of the classical value [Q. There have been attempts to study the effect of electron-electron 
interaction on the noise. For example, noise spectrum in ID systems was calculated from 
Luttinger liquid approximation |5||| . 

In this paper, we study current fluctuations of ultra-small tunnel junctions made of two 
pieces of metal separated by a thin insulating barrier. 

It is helpful to discuss a few limiting cases. It is well known that conventional tunnel 
junctions show Ohmic behavior I(V) = V/Rt when it is driven by external voltage bias and 
the current noise spectrum of Ohmic junctions follows Johnson-Nyquist formula 

^^n^±eV)^(^f). (!) 

Recent advances in fabrication technology have made it possible to construct very small 
devices. It was pointed that electron-electron interaction effect becomes important for small 
devices. For example, it was argued that a tunnel junction develops tunneling gap in current- 
voltage characteristic when the temperature is lower than the charging energy of a junction 
E c = e 2 /2C MM- 

The charging energy was predicted to affect noise spectrum as well. Ben- Jacob et al 
13| studied fluctuations of current in an open tunnel junction. They calculated equilibrium 



noise of a junction and found that the low frequency noise Tiu < Ec is suppressed to zero 
in the low temperature limit k^T <C Ec- 

Following theoretical predictions on the charging effect in a small tunnel junction, many 
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experimental attempts were made to observe the charging effect or so called Coulomb block- 
ade effect. Coulomb blockade effect was verified for mult i-j unction systems |T^-pO|. However 
the clear verification of the effect for a single tunnel junction has not been reported. 

It was later realized that the electrical transport properties of a single tunnel junction 
crucially depends on the electromagnetic environment of a junction which is decided by 
properties of leads on junction chip [pl]-|23l|. Current- volt age characteristic was calculated 
with explicit account of the environment by modeling the environment as linear impedance 
presented to the junction. It was shown that the Coulomb blockade effect is erased by 
the quantum fluctuations of charge when the environment has low impedance and that high 
impedance environment is essential ingredient for the realization of Coulomb blockade effect. 



The dependence on the electromagnetic environment was experimentally observed [p4] 

Behaviors of single tunnel junctions are decided not only by the ratio of temperature 
and the charging energy but also by the environment. It is a goal of this paper to extend 
the earlier works to include noise spectrum calculation for general environment. 

More important goal of this paper is to show that noise measurement in a single tun- 
nel junction is quite nontrivial. For illustration, we calculate current fluctuations at the 
tunneling barrier and fluctuations at the leads. 

Most theoretical works on noise calculated the current fluctuations at a specific point of 
a circuit, usually at a potential barrier right in the middle of a junction. However as pointed 



out by Landauer et al |25[], conventional noise measurements have no relation with the 



electrons crossing the potential barrier of the junction. The noise is measured, for example, 
through the voltage across a series resistor in a circuit. 

Noise measured in experiments is not precisely the same quantity as the noise treated 
in various theoretical calculations. Noise measured in experiments is complicated by the 
intrinsic noise of noise measuring device. To make comparison with theoretical calculations, 
it is usually assumed that the intrinsic noise makes linear contribution and that experiments 
can achieve agreement with theoretical calculations by subtracting off the intrinsic noise 



from the measured noise. We will call this linear superposition assumption [Pq ,p6 



In this paper, we assume that leads, which is connected to a junction and provides electro- 
magnetic environment, can be used as a measuring device and calculate current fluctuations 
at leads. We also calculate current fluctuations at tunneling barrier and compare them with 
current fluctuations at leads to test the linear superposition assumption. To briefly mention 
a result of the calculation, we find that the linear superposition assumption holds only when 



where Z[uS) is the impedance of the environment. In low impedance environment, the 
assumption breaks down due to finite relaxation time of environment. In high impedance 
environment, noise becomes quite nontrivial and spectral density at Tioj « rkEc can be 
reduced by applying small voltage bias. We will later provide the physical origin of this 
novel phenomena. 

In Sec.[II], we present a formulation of a problem which incorporates the charging energy 
and the environment. The formulation allows simple perturbation calculation in the weak 
tunneling limit. In Sec.[TTT], we introduce two current operators. One corresponds to current 
at potential barrier of a junction, which we call tunneling current. The other corresponds to 
current measured at the leads, which we call relaxed current. Expectation values of the two 
operators are evaluated for a few examples of environment. In Sec.[lV|, the noise spectrum 
of the tunneling current is calculated. It shows agreement with existing theories in various 
limiting cases. In Sec.[V|, the noise spectrum of the relaxed current is calculated. It is 
compared to the noise spectrum of the tunneling current and the difference is identified. 
The difference is attributed to the finite relaxation time of environment and nontrivial zero 
point fluctuations of environmental degree of freedom. 



junction is represented by a capacitor notation with an arrow right beside it to emphasize 




(2) 



II. FORMULATION 



Fig.^(a) presents a schematic diagram of a system which we consider in the paper. A 
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tunneling. The junction is connected to the electromagnetic environment and voltage source. 
It is assumed that the environment is linear so that it can be modeled as an impedance Z[lS). 
The impedance is determined by properties of leads connected to junction chip. For detailed 
discussion of environment and the model system in Fig.|l](a), see Ref. ||27|| . 

We use the same formulation as the one adopted in Ref. P%| . The tunnel junction 
connected to the environment is described by the following hamiltonian(for more details 
about this hamiltonian, see P§|), 

H = H qp + H env + H T . (3) 

The first term H qp describes quasi-particle degrees of freedom, 

Hqp = H( e fc + eV)cl a c ka + £ q c\ a c qa , (4) 
ka go- 
where k and q label the quasi-particle states in the left and the right electrodes, respectively 
and a represents spin degree of freedom. 

The second term H env describes environmental degrees of freedom and the charging 
energy, 



\2 N 



(5) 



2C n \e J 2L 

where Q is the charge on the capacitor plates of the junction with the steady state con- 
tribution subtracted off Q = Q — CV . (p is the phase variable across the junction, 
<p(t) = (e/h) jt^ dt'Q(t')/C. Q and (p are conjugate variables describing macroscopic state 
of the junction and they satisfy the quantization relation, 

<p,Q =ie. (6) 

Phase (p is coupled to a set of harmonic oscillators following the Caldeira-Leggett description 
of the environment. 

The third term describes the tunneling between two plates of the capacitor, 

Ht = E T kq c\ a c k(J e-^ + H.c. , (7) 

kqa 



where Tf~ q is the tunneling matrix element. It is simple to verify that e ^ is a charge lowering 
operator, 

e ivQ e -i9 = Q- e . (8) 

Because e ±l ^ explicitly takes care of the charge increase and decrease, the quasi-particle 
degrees of freedom and the environment degrees of freedom can be separated so that quasi- 
particle operators commute with the charge and phase operators. 



III. TUNNELING CURRENT AND RELAXED CURRENT 

To calculate noise spectra, we need to construct current operators. The tunneling current 
operator can be easily identified from the operator equation of motion, 



H, e ^ C \ a C qa 
qa 

■ {iTkqc\ a c ka e~ tip + H.c. 

kqa 



(9) 



This operator represents the net number of electrons tunneled from one plate of the capacitor 
to the other. Experimentally, It corresponds to current measured right at the tunneling 
barrier(Fig.|I|(b)). 

Eq. @ is not a unique choice as a current operator and there is another way of construct- 
ing a current operator, 



~ % r ~ 

Ir= j: \H env , Q 



-Ey-^-Vn)- 



(10) 



For definiteness, we will call It tunneling current and Ir relaxed current. These currents 
can be related in a simple way, 



'R — J-t 



h 



H,Q 



dQ 
dt ' 



(11) 



and from this relation, we interpret In as a charge flow at the environment. Fig.|T](b) illus- 
trates the difference between It and Ir. The construction of the relaxed current operator 
is possible because the environment is explicitly taken into account. 
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Because we are interested in the weak tunneling limit, we take H qp + H env as unperturbed 
hamiltonian and as perturbation. Both H qp and H env are quadratic and the perturbation 
calculation is straightforward. 

Cur rent- volt age characteristic can be obtained from the perturbation theory and one 
obtains 

h{V) = I R {V) = % Im (X ret {-eV/h)) , (12) 
n 

where 

X ret (u) = -t r dtB(t)e lujt (\A(t)e- im ,A\0)e im ]) , (13) 



and A(t) = J2kqa Tk q c \ a (t)ck a (t) ■ (• • -)o represents the average over unperturbed state and all 
operators are written in interaction picture. Because the expectation values of the tunneling 
current and the relaxed current are the same, we will drop the subscript. By explicitly 
evaluating the average, one obtains 

'("') = " ^ L dE T^m Pt - eV ~ • < 14) 

where 

1 r°° 

P ^ = / dt ex P - ^(°))o + iEt / % \ ■ ( 15 ) 

ZTCll J —oo 



Eq.flHJ) together with Eq.(|15|) agrees with the result of Ref. p2| , which justifies two current 
operators. 

All information of the environment is stored in the phase-phase correlation function 
J {t) = ([<p(t) — <p(0)](p(0))o and from the fluctuation-dissipation theorem, J(t) can be related 
to the impedance Z(uj), 

J(t)=2 -^^{coth-— [cos cut- 1]- ism cut \ , (16) 

JO LO Rq [2 J 

where Z t {uj) = {Z~ 1 {uj) — icoC) 1 and Rq is the resistance quantum h/e 2 . 

Let us present current- volt age characteristics in a few special cases. It illustrates the 
effect of the environment and we later find that Si T (u) is closely related to the characteristic. 



Firstly, let us take Z(oj) = 0. Then from Eg. (|15|, [16]), one finds P(E) = 6(E) and the 
junction shows Ohmic behavior. The charging energy plays no role in the low impedance 
environment. Secondly, let us take Z(u) = R ^> Rq and f3Ec ^> 1. Then from Eq. ( |15| , p^ ) , 
one finds P(E) = S(E — Eq) and the characteristic develops a gap and current does not 
flow for \eV\ < Eq- This agrees with the prediction of Likharev et al ||TT| and one realizes 



that their semiclassical theory becomes valid if a junction is connected to high impedance 
environment. Finally, let us take Ohmic environment of arbitrary resistance, Z(u) = R 
30fl . In this case, the characteristic is rather complicated and we present only the limiting 



behaviors at zero temperature. With small voltage bias \eV\ <C Ec, one finds 

\R\eV\ 



1{y) _ exp(-2 7j R/i? Q ) V 



Ro E, 



c 



2R/Rq 

(17) 



T(2 + 2R/R Q ) R T 

where 7 = 0.577... is the Euler constant. Note that Eq.(|17D smoothly interpolates the 
characteristics of the low and the high impedance environments. Also note that with the 
finite impedance, the suppression due to Coulomb blockade effect is not strong enough to 
develop well defined gap even at zero temperature. With large voltage bias eV 3> E c , one 
finds |2Sj the characteristic approaches the high impedance result, 

Rq E 2 c 



i(v) = 4- 

eRi 



eV-E c + 



(18) 



7r 2 ReV 

The zero temperature current-voltage characteristics are presented in Fig.[|. For the exper 
imental measurement of the current-voltage characteristic, see Ref. [|4 



IV. NOISE SPECTRUM OF TUNNELING CURRENT 

A. Noise spectrum of tunneling current 

Noise spectral density of current or noise spectrum is defined as Fourier transform of a 
current-current correlation function, 

s z (u) = \J dtj* Q ({W(o)}) - (i(t))(im) ■ (19) 



Because we later find the noise spectrum of the tunneling current different from the one of 
the relaxed current, we will put subscript T and R to give distinction. 

Calculation of Si t {uj) is as straightforward as the calculation of the current- voltage char- 
acteristic and one finds 



P(hu ±eV-E)(l + e -^ ±eV ^ . 



x 

Structure of Eq.(|20|) is similar to the formula for the current- volt age characteristic Eq.(|14D 
and indeed, Sj t (uj) can be written in terms of I(V), 

Si T M = f E coth g(^±M j (y ± hu/e) . (21) 

Eq.(pTl) allows one to obtain the noise spectrum of the tunneling current from the current- 
voltage characteristic. 

It is a simple task to verify that Eq.([H]) satisfies the fluctuation-dissipation theorem in 
equilibrium V — 0. To verify it, one introduces a coupling of the tunneling current to its 
conjugate force frit), 

H fT = -i T f T (t) , (22) 

and calculates the response function Xt(u) = {Irit))^/ friui). Then, the fluctuation- 
dissipation theorem relates the response to the noise spectrum 

Si T {uo) = h coth ^ Im Xt(^) • (23) 

In the following, we consider equilibrium and excess noise for Ohmic environment at 
zero temperature. Finite temperature noise will be discussed only when it allows a simple 
analytic expression. 

B. Low impedance environment 

Let us first study the low impedance limit Z{uj) = R <C Rq, where the junction becomes 
Ohmic. In equilibrium, Eq. (|2l"D reduces to Johnson-Nyquist noise |yj (FigFJ) 
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S/ T M = coth -— - . (24) 

ttRt 2 

Excess noise is the difference between equilibrium and non-equilibrium noise Sf (u)] V) = 
Si t (lj; V) — Sj t {uj; V = 0). The superscript X is employed to denote the excess noise. 

At zero temperature, the excess noise is absent for \Koj\ > \eV\ and for \hu\ < \eV\, it 
shows piecewise linear dependence on frequency (Fig.|(a)), 

q x \ eV \ ~ IM ^ 

S ^ V )= TrRj* (25) 

Eq.([M|) and Eq.(03j) agree with Eq.([I|) where the charging effect was neglected. The agree- 
ment with the previous study illustrates that the charging effect does not show up in the 
low impedance environment. 



C. High impedance environment 



With Z{u) = > Rq, Eq.(gTj) reduces to 

1 



S It (lu) 
in equilibrium where 



1 + e 



dE 



E 



-f3E 



P(huo - E) 



(26) 



P(E) 



1 



: exp 



(E — E, 



c) 



4:E c k B T 



(27) 



V^E c k B T 

At high temperature k B T > E c , S It {uj) becomes (2k B T/irR T )(l-E c /3k B T+0(E c /k B T) 2 ) 
for Uuj -C Ec, which is the Coulomb blockade correction to white noise. At low temperature 
k B T <C Ec, Si T (uS) becomes {2Ec/'kRt) exp (—Ec/k B T) for hu <C k B T. The Coulomb 
blockade strongly suppresses the noise far below the white noise level in the low temperature 
(Fig-!). 

Both high and low temperature limiting behaviors agree with earlier calculation of Ben- 



Jacob et al where they calculated the equilibrium noise spectrum of an open junction 
which is not connected to the external circuit. Ben- Jacob et a/'s result can be cast into the 
form of Eq.(p6|) with P(E) replaced by 
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Pbj(E) 



ex P 

integer q 



q 2 E c 



£ 6 (E - £ c (2g + 1)) exp ( J . (28) 



Note that Pbj{E) is a discrete form of P{E). P{E) and Pbj(E) become identical in the high 
temperature limit fc^T ^> Eq whereas they exhibit delicate difference in the low temperature 
limit ksT <C Eq. Recalling that P(E) is a probability to emit energy E to environment 
||28|| , one finds that Eq. (|28D implies the quantization of the emitted energy whereas Eq . (f2~7D 
contains continuous spectrum of the emitted energy. 

The charge in each electrode is quantized both in the open junction and in the closed 
junction with high impedance environment. Hence the energy used to charge the electrodes 
has only discrete values. In the open junction, the whole emitted energy is used to charge 
the electrodes. In the junction closed by the environment, however, the emitted energy is 
used only in part to charge the electrodes and the rest excites the environmental degrees of 
freedom. Therefore P(E) has continuous spectrum of energy whereas Pbj{E) allows only 
discrete energy values. 

In the high impedance environment, the excess noise shows piecewise linear dependence 
on frequency at zero temperature(Fig||(a,b)). For small voltage bias \eV\ < E c , the excess 
noise has finite spectral density(Fig.|](a)), 

Z7TKt 

for — \eV\ < \hu\ — Ec < |eV| . Note that even though current does not flow when \eV\ < E c , 
the excess noise does not vanish around hu = ±Ec- At higher voltage bias, finite current 
flows and Fig.|](b) shows the evolution of the excess noise with the voltage. The excess noise 
again shows the piecewise linear dependence. The piecewise linear dependence in both the 
low and the high impedance environments originate from the narrow energy range of P{E). 



D. Ohmic impedance environment 

One can also obtain Si T (u) in arbitrary Ohmic impedance environment Z(u) = R. In 
the low frequency regime \hu\ Eq, the zero temperature equilibrium fluctuations are 
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2R/Rq 



(30) 



suppressed by the Coulomb blockade effect, 

_ exp{-2^R/R Q ) \hu\ InR \hu\ 
It[U) nT(2 + 2R/R Q ) R T [R Q E c 

However the suppression is not complete even at zero temperature. The degree of the 
suppression depends on the impedance and the higher the impedance, the stronger the 
suppression is. 

In the high frequency regime \hu\ 3> Ec, the equilibrium noise eventually converges to 
the result of the high impedance environment, 



Si T (u) 



ttR 1 



M -e c + Rq El 



(31) 



TT 2 R\hu\ 

and the convergence is faster with the higher impedance (Fig||). 

With finite Ohmic impedance, the excess noise no longer exhibits piecewise linear de- 
pendence on the frequency because P(E) allows broad range of energy. Since the Coulomb 
blockade effect is blurred by the quantum fluctuations of charge in the finite impedance en- 
vironment, the excess noise is finite in broad frequency range (Fig.[|(a)). At zero frequency, 



irR\eV\ 

7T T(2 + 2R/R Q ) R t [Rq[~Ec~_ 



sf{u = Q) _leM-2,R/R Q )\eV\ 



1 2R/Rq 



(32) 



E. Noise power 

Before we close the discussion of the tunneling current noise spectrum, let us mention 
one implication of Eq.fl21|). In the presence of finite current, the correlation between charge 
tunneling events can be estimated from the noise power P, which is defined as P = 2ttSi(uj = 
0). If there is no correlation between the tunneling events, the noise power retains its full 
value P = 2el . If there is any kind of correlation which regulates tunneling, it decrease the 
noise power below its full value. For example, in a single channel quantum point contact, 
the Fermi correlation is known to reduce the noise power to P = 2el(l — T) where T is the 
transmission coefficient [0-0]. 
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From Eg. (|2l|) , one finds the noise power of the tunneling current to have the full value 
independent of the environment, 

P T = 2eI(V) , (33) 

where the subscript T is employed to denote that this is the noise power of the tunneling 
current. The full noise power implies that Fermi statistics and the charging energy do not 
give rise to any correlation between tunneling events. This result is understandable since 
we are studying the weak tunneling limit and the time interval between two consecutive 
tunneling events is assumed to be large. 



V. NOISE SPECTRUM OF RELAXED CURRENT 

A. Noise spectrum of relaxed current 

Most theoretical works on current fluctuations calculates the noise spectrum from elec- 
tron crossing of a particular plane in a circuit. However, as Landauer and Martin state K5 



'Currents are measured either through the voltage developed across a small series resistor in 
the circuit, or through the magnetic fields generated by the displacement current. Contrary 
to occasional perceptions, noise measurements have no relation to the crossing, by carriers, 
of a particular plane along the circuit.' 

To clarify the relation between conventional noise measurements and noise spectrum 
treated in various theoretical works p[,f7 ,pL3 , 32|| , it is necessary to take an explicit account 



of a measuring device in the formulation. We assume that leads connected to a junction 
can be used as a measuring device and calculate current fluctuations at leads. The relaxed 
current Ir corresponds to charge flow at leads and therefore Si r (uj) corresponds to spectral 
density of current fluctuations measured at leads. We will study the difference between 
Si R (u) and Si T (u) in this section. 

No tunneling limit Rt — > oo illustrates a crucial difference of Si R (u) and Si T (u). In 
this limit, the tunneling current operator itself vanishes because it is proportional to the 
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tunneling matrix elements. Hence, one finds Si T (u) = 0. This result is natural since It is 
the electron flow measured right at the tunneling barrier and electrons cannot tunnel the 
barrier at all in the no tunneling limit. On the other hand, the relaxed current operator 
does not depend on the tunneling matrix elements explicitly and in the zero tunneling limit, 
one finds Sj r (uj) = Sj° r \uj) where 

S<® H = -Re ( * - ) hu coth ^ . (34) 



7T 



Z(u) - 



The superscript (0) denotes that this is the zeroth order contribution in the perturbation 
calculation. Note that Sj R \u) is a generalization of Johnson-Nyquist noise |3l| Sj(u) = 
\/{itR)Tiuj coth(ptkj/2). Si R (u) measures the fluctuations of the relaxed current in the 
environment and Si R (u>) does not vanish even in the no tunneling limit because of the zero 
point fluctuations in the environment. 

With finite tunneling, it is usually assumed that Si a (u) is a linear superposition of the 



intrinsic noise S\° r (uj) and the tunneling noise Sj t (lu), which we call linear superposit 
assumption. The assumption is frequently used to make comparison of experiments with 
theoretical calculations (for example, see Ref . |^6| ) . 

Explicit calculation of Si R (u) to the second order in the tunneling matrix produces 

sv>) = s£V) + sffH (35) 



ton 



where 



^'^ MWf 5 "'" 1, (36) 

S?^\uj) = — Trim i 77 I coth — — 

xRe Y, {Xret(-eV/h ±u)- X ret (-eV/h)} . 
± 

The superscript (2) in Sj~j(u) denotes that this is the second order contribution in the 
perturbation calculation. Unlike the intrinsic noise, S Ir (lj) depends on the voltage. 
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At zero frequency, one finds Sj^(O) = Si T (0) and Sj R (0) = Sj 2 ^\o) = 0, so that 
that the linear superposition assumption holds at zero frequency for arbitrary environment. 
Because Sf^ (00) vanishes at zero frequency, one also finds that noise powers of the tunneling 
current and the relaxed current are always the same Pr = Pp. For non-zero frequency, one 
finds that the assumption holds if 

|^)| «-^|, (37) 

and it does not otherwise. Eq . (p7|) has interesting implications. For example, in the Ohmic 
environment with high impedance Z(u) = R^> Rq where Coulomb blockade effect is clear, 
the excess noise of the tunneling current exhibits interesting feature at frequencies around 
hu) = ±Ec- However Eq. (|37j) predicts that at this frequency range, the excess noise of the 
relaxed current is noticeably different from the tunneling excess noise since Eq. (|57| ) translates 
to \hu\ < (R Q /R)E c (see Sec.[VT]). 

Let us first discuss the equilibrium noise. Due to the fluctuation-dissipation theorem, 
the equilibrium noise spectrum can be alternatively derived by coupling the relaxed current 
to its conjugate force 

H fR = -I R f R (t) , (38) 



and calculating the response function. This procedure produces 

?< >^-jl coth^ 

7ch 2 2 { (1 - iuCZ (w)) 

x (x ret (u) + x* ret {-u) + -x ret (o) - x; et (o))} . 



SiM = Sf'ico) - coth ^Im { - -^rr^ (39) 



It can be verified that Eq. (|39|) agrees with Eq.(|35|) in equilibrium. Note that though both 
Si R {ui) and Si T {u) satisfy the fluctuation-dissipation theorem, they are still different since 
they are related to different response functions. 



B. Low impedance environment 



To make predictions from Eq.(^), one has to evaluate the real part of the retarded 
Green's function X ret (uj). For uniform tunneling density of states, one finds at zero temper- 
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ature 



Re {X ret (-eV/h ±lu)- X ret (-eV/h)} (40) 
± 

h 2 



JdE(P(E)-P(-E)) 



2ire 2 R T 

x Yl H eV -E±hu) In \eV — E ± Tiu\ - (eV - E) In \eV - E\} . 
± 

In the low impedance environment Z(uS) = R Rq, one finds S\ 2C \uj) vanishes at zero 
temperature and the zero temperature equilibrium noise becomes 

S lR {u) = + -^Re ( , 1 ) \huj\ . (41) 

ttR t \ (1 — iuCZ(u)) ) 

It can be verified that Eg. (flip agrees, to the leading order in 1/Rt, with zero temperature 
equilibrium noise in the presence of a shunt resistor Rt in parallel with a capacitor C 

sm= > U) + w ) ih - (42) 

When the voltage bias is applied, one finds the zero temperature excess noise 

s ^) = T^mr^M- (43) 

Eq.(^3]) can be understood if we assume that current fluctuations induced by tunneling is 
incoherent with equilibrium zero point current. 

Then the behavior of a junction is essentially classical. The only effect of the tunneling is 
to introduce additional charge to plates of capacitor and from a conventional circuit analysis 
of a capacitor connected to impedance Z(u) |j33| , one finds the current induced by tunneling 

I R (t) = J dsTZ(t - s)I T (s) , (44) 

where TZ(t) = for t < and 1Z(u>) = 1/(1 — iuCZiuj)). From the incoherence assumption, 
one immediately finds 

s ^ = \i-Jcz M \* SfM ' (45) 
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which agrees with Eq.fl43|). 

The factor 1/|1 — iuCZ(uj)\ 2 in Eq.(fi5|) originates from the finite relaxation time of the 
environment; the relaxation of tunneled charge to the environment is not immediate and it 
takes characteristic time of CZ(uj). 

Fig.^(a) shows the relaxed excess noise in comparison to the tunneling excess noise. At 
zero temperature, there are three relevant energy scales eV,h/RC, Eq- Fig.||(a) assumes 
eV ^> h/RC ^> Ec- Due to the low impedance of the environment, the junction shows 
Ohmic behavior. The tunneling excess noise(dashed line) decays linearly with the frequency 
and its peak width is eV. The relaxed excess noise(solid line) has the same magnitude 
with the tunneling excess noise at zero frequency. However, its peak width is given as 
min(eV, h/RC) and for the case considered in Fig.[|(a), it is h/RC, making the decay of the 
relaxed excess noise faster than the tunneling excess noise. 

C. High impedance environment 

In the high impedance environment, the junction becomes highly nonlinear and the shunt 
resistor model cannot explain its behavior. To understand the junction behavior in the high 
impedance environment, the investigation of the zero point fluctuations in the environment 
is needed. 

In the no tunneling limit, the environment is decoupled from quasi-particle degrees of 
freedom. The environmental hamiltonian is quadratic and therefore it can be diagonalized 
as a set of independent harmonic oscillators. Each harmonic oscillator shows zero point 
fluctuations of the magnitude (x 2 ) = h/ (2m uj ) coth j3hu /2 where m and u are mass and 
frequency of the oscillator, respectively. Harmonic oscillator with frequency ujq contribute to 
the spectral density of the relaxed current at uj = ujq and the zeroth order spectral density 
Sj°J (on) is a sum of such contributions. 

The situation becomes complicated if one turn on tunneling. The tunneling couples 
each harmonic oscillator to quasi-particle degrees of freedom and other harmonic oscillators. 
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Therefore the tunneling provides each harmonic oscillator with the coupling to an effective 
heat reservoir which consists of other harmonic oscillators and quasi-particle degrees of 
freedom. 

The effect of such coupling is twofold. It provides a noise source and it also introduces 
dissipation. The noise source tends to increase fluctuations in, for example, harmonic oscil- 
lator displacement variables while the dissipation has the opposite effect. It is known (for 
example, see Ref. ||34| ) that the coupling to the reservoir can even reduce the zero point 
point fluctuations in certain situations. 

The relaxed current Ir is linear in the environmental variables. Therefore it can be 
written as a linear combination of harmonic oscillator displacements and the fluctuations of 
Ir is directly related to the fluctuations of harmonic oscillators. 

A harmonic oscillator coupled to the reservoir is a difficult problem in general. However 
in the weak tunneling limit, one can perturbatively calculate the effect of the reservoir. 

It is worth mentioning that the voltage bias is a variable describing quasi-particle degrees 
of freedom. Therefore fluctuations of harmonic oscillators and the relaxed current depend 
on the voltage bias. 

Fig.^(b) shows the relaxed excess noise in comparison to the tunneling excess noise at 
zero temperature in the high impedance Ohmic environment Z(u) = R ^> Rq. The relaxed 
excess noise can be interpreted as change of harmonic oscillator fluctuations by the voltage 
bias. It is assumed that E c > eV ^> h/RC. Note that current does not flow because 
the voltage bias is below the gap eV < Ec- The tunneling excess noise shows two peaks 
at hu = ±Eq- In comparison, the relaxed excess noise (solid line) becomes negative at 
hu ~ ±Ec- The relaxed excess noise is multiplied by 10 4 to magnify the negative dips. 
Though the dips are quite shallow, it is important to notice that the dips imply the decrease 
of the total noise upon the application of the voltage bias. The dash-dotted line shows 
1/(1 + u 2 R 2 C 2 ) S? R (uj). The dash-dotted line is also multiplied by 10 4 . At higher voltage 
bias eV > Ec, the relaxed excess noise shows a peak at hu = and shallow dips at 
hu ps ±E C - 

18 



From Fig.||(b) and the relation Sf(0) = Sf T (0), one finds that the dissipative effect is 
dominant at Tiuj « ±E C while the noise source effect is dominant at fhu ~ 0. 



D. Ohmic impedance environment 

Fig.^(c) shows the two excess noises at zero temperature in the Ohmic environment 
Z(uj) = R. It is assumed that R = Rq and eV = 0.5Ec- The tunneling excess noise(dashed 
line) shows two peaks, while the relaxed excess noise(solid line) shows a single peak at 
hu = 0. The dash-dotted line shows 1/(1 + uj 2 R 2 C 2 )Sf T {uj) for comparison. The dash- 
dotted line exhibit qualitatively the same feature as Si R (u) though there are quantitative 
difference. 



E. Experimental implications 

Eq.(|35|) makes one interesting prediction in the low impedance Ohmic environment. As 
illustrated in Sec. |V B| (Fig.[5](a)), the excess noise measured from the environment is different 



from the tunneling excess noise and follows Eq.( 

Recently, there was an experiment which measured the excess noise of the quantum point 
contact in the low impedance environment ||35|| . The voltage bias V of the order of lmV was 
applied and the excess noise was measured at frequency uj ~ 10GHz which is much lower 
than the voltage hu <C eV. The frequency was chosen to avoid complications due to other 
noise sources at low frequency such as 1/f noise. The measured excess noise was compared 
with theoretical calculations of the noise power (for example, Ref. |32| or S^,(u;)) based on 



the assumption that the non-zero frequency excess noise is essentially the same as the noise 
power if Tiuj <C eV. The experiment showed nice agreement with existing theories in ballistic 
regime. However in the pinched off regime which corresponds to the weak tunneling limit, 
it was found that the noise was reduced to about 1/3 of theoretical predictions. 



The quantum point contact in Ref. [53 is induced electrostatistically in the plane of 



a 2DEG embedded in GaAs-AlGaAs heterostructure. Due to the very different density of 
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charge carriers in metals and semiconductors, our calculation on metallic junctions cannot be 
directly applied to analyze the experiment. However we believe that the difference of Sf (u) 
and Sf(u) in the low impedance environment is a generic property of the weak tunneling 
limit because the difference originates from the finite relaxation time of the environment 
which does not depend on charge carrier density. 



We argue that the noise suppression observed in Ref. |35[ may have its origin in the 
different frequency dependence of the relaxed and tunneling excess noise as illustrated in 
Sec. [V B| . However, we are unable to establish quantitative agreement of the experiment with 
Eq.(f43D because the precise values of R and C are not available. More experiment on, for 
example, the frequency dependence of the suppression factor, is needed. 

Eq.(|35|) makes another interesting prediction in the high impedance Ohmic environment. 
At frequency hu « ±E C , the spectral density Si R (u) measured at the leads decreases below 
its equilibrium value when the voltage bias is applied ( Sec , |V Cp , making the relaxed excess 
noise negative. Usually the spectral density increases when the voltage bias is applied. 
Therefore the observation of the negative excess noise can serve as a crucial test of the 
present calculation. 



F. Weak scattering limit 

Calculation in this paper is restricted to the weak tunneling limit Rt ^> Rq- It will be 
interesting to carry out similar study for the weak scattering limit. According to a recent 
experiment on quantum point contact []35|| , the measured excess noises show the deviation 
from existing theories only in the pinched off regime and as a channel opens(as a transmission 
coefficient T approaches 1), the deviation disappears. If the deviation is due to the difference 
between the relaxed excess noise and the tunneling excess noise, the experiment implies that 
the junction in the weak scattering limit no longer behaves as a conventional capacitor and 
the difference between the two excess noise vanishes. Though the experiment used junctions 
in semiconducting heterostructure, we expect metallic junctions show similar behavior at 
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least in the low impedance environment. 

Also there was a theoretical study of a single tunnel junction using the functional integral 
technique ||36|| . It was claimed that the junction becomes purely Ohmic and the capacitive 
effect vanishes in the weak scattering limit. 

In the weak scattering limit, electron transport is continuous rather than being sporadic 
as in the weak tunneling limit and we don't expect a conventional capacitor is a proper 
description of the junction. However, further study is required to gain the quantitative 
understanding in the weak scattering limit. 



VI. CONCLUSION 

In summary, we considered current transport through a junction connected to electro- 
magnetic environment in the weak tunneling limit. The charging energy is taken into account 
to incorporate the electron-electron interaction effect. We calculated both the tunneling 
noise and the relaxed noise for arbitrary environment and showed their dependence on the 
charging energy and environment. We found that contrary to conventional expectation, the 
current through the potential barrier of the junction and the current measured from the 
environment have different excess noise. In the low impedance environment, the difference 
is due to the finite relaxation time of the junction. In the high impedance environment, 
additional complication occurs and at frequency Tiu « ±Ec, the spectral density of the 
relaxed current decreases below the equilibrium level when small voltage bias is applied. 
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FIGURES 

FIG. 1. (a) Schematic diagram of the circuit. The circuit contains a single tunnel junction 
, voltage source, and the impedance Z(oS) of the environment, (b) The tunneling current and 
the relaxed current. The tunneling current represents the electron flow right at the potential 
barrier of a junction. As an example, the tunneling current through the insulating barrier in 
metal- insulator-metal junctions is drawn. The relaxed current represents the electron flow through 
the environment which is connected to a junction in series. In the figure, Ohmic environment with 
resistance R is drawn. 

FIG. 2. Zero temperature current- voltage characteristics with Z(uS) = 0.01i?Q (solid line), 
Z(lo) = WORq (dashed line), and Z(u>) = Rq (dash-dotted line). 

FIG. 3. Zero temperature equilibrium noise spectrum of the tunneling current for 
Z(lo) = O.OIRq (solid line), Z{uj) = WORq (dashed line), and Z{uS) = Rq (dash-dotted line). 

FIG. 4. (a) Zero temperature excess noise spectra of the tunneling current for 
Z{uS) = 0.01i? Q (solid line), Z(cj) = 100i?Q(dashed line) and Z(uj) = Rq (dash-dotted line). Volt- 
age bias of eV = 0.5Ec is assumed, (b) Evolution of zero temperature tunneling excess noise 
with voltage bias in the high impedance environment Z{uS) = WORq. eV = 0.5i?c(solid line), 
eV = l.5Ec (dashed line), and eV = 2. 5Ec (dash-dotted line). 

FIG. 5. The excess noise spectrum of the relaxed current Sf R (w) (solid line) vs. the excess noise 
spectrum of the tunneling current Sf T (uo) (dashed line) at zero temperature. The dash-dotted line 
shows 1/(1 +uj 2 R 2 C 2 )Sf T {uj) for comparison, (a) Z{w) = OMRq, eV = 200E C . The dash-dotted 
line is not visible because it overlaps with the solid line, (b) Z{oj) = 100-Rq, eV = 0.5Ec- The 
solid line and the dash-dotted line are multiplied by 10 4 to magnify their features. Note that 
Sf R {u) is negative, (c) Z(u>) = Rq, eV = 0.5Eq- Sf R {uo) has a single peak at u> = compared to 
two peaks of Sf T {uj). Sf R (uj) deviates from 1/(1 + uj 2 R 2 C 2 )Sf T {uj) due to nontrivial behavior of 
the zero point fluctuations. 
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